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We introduce a path to a possibly new class of magnetic materials whose properties are determined 
entirely by the presence of a low concentration of specific point defects in their crystal structure. 
Using model Hamiltonian and ab-initio band structure methods we demonstrate that even large 
band gap nonmagnetic materials as simple as CaO can exhibit extraordinary properties like half 
metallic ferromagnetism upon introducing a small concentration of Ca vacancies. We show that 
such defects will initially bind the introduced charge carriers at neighboring sites and depending on 
the internal symmetry of the clusters formed by neighboring sites form "local" magnetic moments 
which for concentrations as low as 3% transform this non-magnetic insulator into a half metallic 
ferromagnet. 



The recent discovery of ferromagnetism with high tran- 
sition temperatures and very small magnetic moments 
in the hexaboride s |l| , Co substituted T1O2 and Co 
substituted ZnO jjjhas opened a lively discussion re- 
garding the role of defects in transforming insulating 
non magnetic compounds into ferromagnets with inter- 
esting properties In the hexaborides Monnier and Delley 
demonstrated, using density functional band structure 
approaches, that neutral Bg vacancy's in a superlatticc 
can result in a ferromagnetic ground state even for a low 
density of such defects We note that the nominal 
valence of the clusters in fairly ionic materials such 
as CaB 6 is 2— and so two electrons are needed to charge 
compensate for such a B6 vacancy. It is also interesting 
to note that substituting divalent Co for a tetravalent Ti 
in Ti02 also can be compensated with two holes in the 
Oxygen 2p valence bands which could be bound impurity 
states in close proximity to the Co impurity. These obser- 
vations may remind some that in hydrocarbon molecules 
with ring structures Longuet-Higgins predicted magnetic 
ground states due to a kind of molecular Hund's rule 
coupling well known in atoms |!|. This property has in 
fact been used in the past by Torrance et al Q in at- 
tempts to make organic molecular ferromagnets. More 
recently Eskes has used a Hubbard model calculation of 
ring systems and clusters to demonstrate again a kind 
of molecular Hund's rule coupling for orbitally degener- 
ate states, leading to very stable magnetic (spin triplet) 
ground states for either two electrons or two holes, in 
particular geometries 0. 

In the present work we report on the influence of di- 
lute divalent cation vacancies in oxides with the rock salt 
structure. We predict that such systems have ferromag- 
netic ground states with small magnetic moments but 
possibly high transition temperatures. For 3.125% Ca 
vacancies in CaO LDA band structure calculations pre- 
dict a half metallic ferromagnetic ground state which if 



correct could play an important role in the modern field 
of spintronics and opens a path to a whole new class 
of ferromagnetic materials. Before presenting the band 
structure calculations we discuss the physics we believe to 
be responsible for the magnetic behaviour found. In the 
rock salt structure the cations are surrounded by an octa- 
hedron of anions with filled valence p bands as shown in 
Fig.|l|. As also shown, the a-bonding anion p orbitals with 
lobes directed to the cation site, would be the most en- 
ergetically favourable orbitals to accommodate the holes 
needed to charge compensate a cation vacancy at the 
center of such an octahedron. Because of the negative 
effective charge of the cation vacancy the holes in the an- 
ion valence p band feel a strong attraction to the vacancy 
site and if large enough this attraction can bind two holes 
as in fact is found in the LDA calculations. 

In order to get insight into the electronic structure of 
such an octaheral cluster of 2 ~ ions each with a full 2p 
shell we use a tight binding like molecular orbital (MO) 
approach. The symmetries of the MO's formed by the 
cr-bonding orbitals are di g , t\ u , and e g at energies equal 
to — Atpp + e a , e a , 2t pp + e a respectively. Note that simple 
sum of orbitals shown in Figj^(a) represents the a\ g sym- 
metry state and this state has the lowest energy for elec- 
trons and the highest for holes. The 7r-bonding orbitals 
form the symmetries t2 g , fi«, t2 U , and t\ g and have en- 
-2t pp7r , and 2t pp respectively where 
tppix) with the a and ir hopping inte- 
grals now relative to the 0-0 bond direction and given 
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Note that for simplicity the zero of energy is chosen to be 
at e^. One can see that the irreducible representation t\ u 
appears in both the a and tt combinations. As a result 
the molecular orbitals of only this particular symmetry 
will mix with each other resulting in a shift of their ener- 
gies to l/2e CT + t ppT 
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FIG. 1: An artists concept of the oxygen a (a) and tt (b,c) 
bonding orbitals relative to the O-vacancy bond direction sur- 
rounding a Ca vacancy. Also shown are the definitions of the 
hopping integrals (b) t pp , (c) t pp7T , and (d) t' pp given in terms 
of electrons. 



The hopping integral t' pp is defined as l/2(£ pp(T +tp Fn .) and 
is shown in Fig. 0(d). 

The schematic energy diagram is given in Fig. | The 
lowest energy states for holes are the e g states formed 
by a linear combination of the p orbitals with lobes di- 
rected to the vacancy with phases such as to result in a 
e g symmetry state when viewed from the vacancy and 
the tig states with a corresponding linear combination 
of p orbitals with lobes directed perpendicular to the 
O- Vacancy bond direction. A simple calculation shows 
that a hole in p orbitals with lobes pointing towards the 
vacancy has an energy about leV lower than a hole in 
a p-orbital with lobes perpendicular to the O- Vacancy 
bond direction because of the crystal field produced by 
the effective 2— charge of the cation vacancy. There- 
fore the lowest energy single hole state will be the doubly 
degenerate e g molecular orbital (Fig. ^(a)). We note 
that for anion vacancies the conduction band orbitals of 
the cation would be occupied by the charge compensat- 
ing electrons and their energy would be lowest for the 
non-degenerate a\ g orbital as also shown in Fig.||(b). 

Having understood the single-particle molecular or- 
bitals, we move to the case of two holes (electrons) in 
the completely filled (empty) shell formed by the anion 
p (cation s) orbitals in the case of cation (anion) vacan- 
cies. If the Coulomb interaction between those particles 
is negligible, then two holes will occupy the e g molecular 
orbital and they will form three degenerate states 1 Ai g , 




FIG. 2: Schematic single particle energy level diagram for 
(a) HOLES in anion orbitals and (b) ELECTRONS in cation 
orbitals. Dash line in fig. (a) represents a shift of energy levels 
of molecular orbitals with t\ u symmetry due to the mixture 
between a and 7r-bonding anion orbitals. The energy of the 
atomic a and n orbitals are indicated by arrows on the left 
hand margine. 

1 E g and 3 ^2g (Fig.||(b)). Two electrons on the other 
hand will result in a non degenerate singlet ground state 
(Fig.||(a)). Let us now switch on the Coulomb repulsion 
of which the dominant term will be that when two par- 
ticles are on the same site as in a Hubbard model. The 
three states found for the configuration with two holes 
will now split up because for the singlet states the two 
holes can be in the same site and in the same a orbital 
while this is forbidden for the spin triplet state. This 
means that the singlets will increase in energy leaving us 
with a triplet and magnetic ground state. From Auger 
spectroscopy studies of a large number of simple oxides 
the on-site cooulomb repulsion energy of two holes in an 
oxygen p orbital is 5-7 eV [|[ || ■ 

To shed more light on this issue we carried out ex- 
act diagonalization calculations for the single octahedral 
cluster with two electrons in s orbitals as well as with 
two holes in the p orbitals of a octahedral cluster taking 
into account the on-site Coulomb interaction U exactly, 
and the results are shown in Fig. [| One can immedi- 
ately see that while a spin singlet state is the ground 
state for two electrons independent of the value of U, a 
spin triplet state is the ground state for the two holes as 
soon as U is different from zero. This triplet state has 
symmetry 3 A2 g and the wave function consists mainly of 
a (e g ) 2 configuration shown schematically in Fig. ||(a). 
Note that the lowest singlet state is at about 0.3-0.4 eV 
higher in energy than the triplet ground state of the two 
holes for realistic parameters. As mentioned above this is 
a kind of molecular Hund's rule coupling and is a result 
of the degenerate nature of the lowest energy molecular 
orbital for holes in a octahedral cluster. We note that 
for oxygen vacancy's it does not matter if we are dealing 
with s or p, d states of the cations. As long as band is 
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FIG. 3: Energy diagrams of the lowest three states for two 
particles in an octahedral cluster as a function of the on-site 
Coulomb repulsion U; (a) Two electrons introduced into the 
empty s-orbitals with the hopping integral t ss = —1.0 eV and 
(b) two holes introduced into the completely filled p-orbitals 



with the hopping integrals t p 



1.2 eV and t v 



-0.3 eV 



and the crystal field potential e a — 1 eV and £„ = eV given 
in terms of electrons. Solid (open) symbols are for spin triplet 
(singlet) states. 



empty the lowest energy state for electrons in octahedral 
geometry will always be a totally symmetric state and 
therefore non degenerate. 

Now that we have found that cation vacancies in sim- 
ple divalent monoxides with the rock salt structure or 
in general cation vacancies in compounds with octahe- 
dral coordination form impurity states with local mag- 
netic moments we may ask what the magnetic state of 
the material is for a finite concentration of such vacan- 
cies. Do they couple ferromagnetically or antiferromag- 
netically and is the resulting material metallic or still 
insulating? To address this we consider the problem of 
a super lattice of cation vacancies with a 2x2x2 super- 
cell. This would correspond to a vacancy concentration 
of only 3.125%. We calculate the electronic structure us- 
ing a density functional band structure code (TBLMTO- 
47 ) based on the local density approximation (LDA) . 
For the Ca atoms, the basis from 4s, 4p, and 3d orbitals 
was used, while for oxygen, the 2s, 2p and 3d states were 
considered. The electronic states inside the empty sphere 
between the atoms were expanded up to 3c? states includ- 
ing Is and 2p orbitals. 

Before we discuss the results of band structure calcu- 
lations we note that, due to the ionic character of this 
compound, the valence band is formed of mainly oxygen 
2p orbitals and is full and the conduction band is formed 
of mainly Ca orbitals and is empty in the vacancy free 
compound. Pure CaO is therefore found to be a diamag- 
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FIG. 4: Total density of states of pure CaO (upper panel) and 
compound with vacancy in anion (second panel) and cation 
(two lower panels) sublattices. The zero of energy is at Fermi 
energy. 



netic insulator with a calculated band gap of 3.46 eV. 
The experimental gap is 7 eV |llj . 

In Fig. [| we display the total density of states of pure 
CaO and CaO with vacancies in the cation or anion sub- 
latticies. We observe the impurity state just below the 
conduction band and the chemical potential shifted up 
above this impurity state for the case of Oxygen vacan- 
cies. This material is a non magnetic metal. The same 
calculation for MgO with Oxygen vacancies exhibits a 
clear gap separating the defect induced band from the 
conduction band resulting in an insulating non magnetic 
ground state. The details of calculations on MgO, CaO, 
and SrO will be presented elsewhere. 

In sharp contrast the band structure for Ca vacancies 
exhibits also well defined defect state bands crossing the 
Fermi energy resulting in a metallic solution. Also here 
we see the clear split off character of the vacancy induced 
states consistent with the physical picture we presented 
above of the effects of the defect potential. Even more 
interesting is the fact that the Ca vacancy material is 
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FIG. 5: Ca vacancy partial density of states. The zero of 
energy is at Fermi energy. 

found to be ferromagnetic and the bands crossing the 
Fermi energy are totally spin polarized. A Ca deficiency 
results in a magnetic moment of 2fj,s per super cell, 88% 
of which is concentrated on the six Oxygen ions that are 
the nearest neighbors to the Ca vacancy. 

In both cases of Ca or oxygen vacancies we find that 
the charge compensating holes or electrons are rather 
strongly bound to the vacancy. In Fig. || and Fig. ^ 
we show the partial density of states projected on the 
impurity site. It is easy to see that in the case of Ca va- 
cancies the holes are mainly d like with e g symmetry but 
band structure effects also introduce quite a strong com- 
ponent of p symmetry at the Fermi energy resulting in a 
half-metallic ferromagnet. On the other hand, electrons 
induced by oxygen vacancies occupy local molecular or- 
bitals of s symmetry and result in a nonmagnetic ground 
state. Both of these agree very well with the results of 
our cluster calculations. We suggest that the picture pre- 
sented describes a path to new classes of magnetic mate- 
rials and is a very general phenomenon. 

The general underlying rule for obtaining a magnetic 
ground state is that the crystal structure must be such 
that the ground state of the charge compensating molecu- 
lar orbital is orbitally degenerate. Secondly the local im- 
purity potential should be large enough to quite strongly 
bind the charge to the nearest neighbor atoms which is 
needed to validate the treatment as a cluster or molecule 
which then will result in the formation of a "local" mo- 
ment. An on-site Coulomb interaction will then result in 
a high spin ground state for more than one electron or 
hole in such a molecular orbital which will mediate a fer- 
romagnetic magnetic coupling between such clusters via 
either a superexchange or a double exchange like mech- 
anism. We note that for the case considered above the 
charge compensating clusters have considerable overlap 
even though the concentration is only a few per cent. We 
also note that such degenerate ground states are common 
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FIG. 6: Density of states projected onto the O vacancy site. 
The zero of energy is at Fermi energy. 



for high symmetry systems such as ring like structures in 
addition to the octahedron considered here. 

In conclusion we have demonstrated that point struc- 
tural deformations of crystals such as vacancies can in- 
deed confine the compensating charges in molecular or- 
bitals formed by atomic orbitals on the nearest neigh- 
bours. We have shown that under certain conditions local 
magnetic moments will be formed due to a kind of molec- 
ular Hunds rule coupling with energetics determined by 
kinetic energy and symmetry considerations rather than 
exchange interactions. This can lead to high temperature 
ferromagnetic ground states and in the example discussed 
even a half metallic ferromagnetic ground state. Because 
of the rather general nature of the physics presented we 
believe that this may point to a path towards a new class 
of ferromagnetic materials. 
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